Nuclear spin ferromagnetic phase transition in an interacting 2D electron gas 
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Electrons in a two-dimensional semiconducting heterostructure interact with nuclear spins via the 
hyperfine interaction. Using a a Kondo lattice formulation of the electron-nuclear spin interaction, 
we show that the nuclear spin system within an interacting two-dimensional electron gas undergoes 
a ferromagnetic phase transition at finite temperatures. We find that electron-electron interactions 
and non-Fermi liquid behavior substantially enhance the nuclear spin Curie temperature into the 
mK range with decreasing electron density. 
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The use of the electron spin as a qubit for quantum 
computing relies on the ability to coherently control sin- 
gle electron spins in semiconductor quantum dots 
Over the last years much progress has been made for dots 
in GaAs semiconductors, where single spin lifetimes have 
been measured to range far into the ms-range @, 0, 0], 
and where coherent manipulation of single- and two-spin 
states was successfully implemented 0, @|. Still, a major 
obstacle to further progress is the comparatively short 
spin decoherence time in these materials, ranging from 
100 ns in bulk [3] to fis in dotsjl|. The main source of 
decoherence for a single electron spin confined to a GaAs 
dot is coming from the contact hyperfine interaction with 
the surrounding nuclear spins @,S Several ways to 
overcome this problem have been proposed such as spin 
echo techniques [H,0], projection of the nuclear spin state 
0] or polarization of the nuclear spins 0, H, [l(J EH] . How- 
ever, in order to extend the spin decay time by one order 
of magnitude, a polarization of above 99% is required 
0, which is still far away from the 60% so far reached 
in quantum dots via optical pumping [12]. One way to 
overcome this problem would be that nuclear spins be- 
come fully polarized at low enough temperatures, with- 
out any external magnetic field or optical pumping. This 
is the case if the nuclear spins undergo a ferromagnetic 
phase transition at a finite Curie temperature T c . Quite 
remarkably, the possibility of such a nuclear-spin phase 
transition to occur in a metal was studied more than sixty 
years ago by Frohlich and Nabarro (FN) Using a 

Weiss mean field argument they showed that the Curie 
temperature T c of nuclear spins in a three dimensional 
(3D) metal becomes 
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where A denotes the hyperfine coupling strength between 
the nuclear and electron spin and Ef the Fermi energy. 
For a typical metal, T c is of the order of micro-Kelvin 
or less. However, for a two-dimensional electron gas 
(2DEG) in GaAs semiconductors, Eq. (TJ would pre- 
dict nuclear ferromagnetism with T c ~ lmK, which is 



surprisingly high 

However, the direct use of Eq. (TJ, which was derived 
for a bulk metal, to a 2DEG in a semiconductor is very 
problematic. The purpose of this letter, therefore, is to 
reconsider this issue for a 2DEG and to estimate the nu- 
clear spin Curie temperature. Our analysis below will be 
based on the Kondo lattice model , where we integrate 
out the electron degrees of freedom to derive an effective 
spin Hamiltonian whose exchange is given in terms of the 
static electronic spin susceptibility Xs(l)- Using a spin- 
wave analysis, we will show that the electron-electron (e- 
e) interactions in the 2DEG and the induced non-Fermi 
liquid behavior in Xs(o) (l5l. fl^. [TtI. fl8| ultimately enables 
a ferromagnetic phase transition of the nuclear spins. 
For sufficiently strong interactions and/or low electronic 
densities (with the dimensionless interaction parameter 
r s ~ 5 — 10) the Curie temperature can be pushed into 
the milli-Kelvin regime, and thus, the phase transition 
should become accessible experimentally. 

Model Hamiltonian. In order to study an interacting 
2DEG coupled to nuclear spins within the 2DEG, we 
adopt a tight-binding representation in which each lat- 
tice site contains a single nuclear spin and electrons can 
hop between neighboring sites. The Hamiltonian describ- 
ing such a system reads 



H = H 



1 



N 



2 ^ 



A 



jC ja v al 3C jl3 ■ I j =H + H n 



(2) 



where H denotes the conduction electron Hamiltonian 
and H n the electron-nuclear spin hyperfine interaction. 
Hq can be rather general and includes e-e interactions. 
In Eq. ([2]), Cj a creates an electron at the lattice site 
fj with spin a and a represent the Pauli matrices. We 
have also introduced Ij the nuclear spin located at the 
lattice site fj, and Aj the hyperfine coupling constants 
between the electron and nuclear spins at site fj. The 



electron spin operator is defined by Sj = 



2°j'a' 



N denotes the total number of sites on the 2D lattice. 
In our formulation, the nuclear spin density is n s = a~ 2 
where a is the lattice spacing. From here on, we assume 
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Aj = A > which means we assume the hyperfine in- 
teraction to be the same for all atoms that constitute 
the heterostructures (typically Ga and As). We also ne- 
glect direct dipolar interactions between the nuclear spins 
which is in general smaller than the indirect interaction 
as we will see. This amounts to assume that the dipo- 
lar interaction energy scale Edi P is among the smallest 
one and particularly that feT > Edi p , where T is the 
temperature. This assumption is crucial since it allows 
us to focus on the nuclear spins which are within the 2D 
electron gas thickness (in growth direction) and justifies 
our 2D description [19]. 

The general Hamiltonian in Eq. |(2]) is the well-known 
2D Kondo lattice Hamiltonian (KLH), though Hq con- 
tains also e-e interactions. The regime we are interested 
in corresponds to the weak Kondo coupling regime in the 
sense that A -C Ep, where Ep is the Fermi energy. The 
KLH has been introduced to describe various physical 
properties of heavy-fermion materials [3, H3], and more 
recently also of ferromagnetic semiconductors [21]. 

Before turning to the extended system let us briefly 
consider the special case of a single electron confined to 
a quantum dot which interacts typically with 10 6 nuclear 
spins 0, [l3]. This case can be described by the above 
KLH by allowing in Hq for a confinement potential for 
the dot, which provides the largest energy scale. Indeed, 
we can then project H n into the ground state of Hq, and 
the hyperfine Hamiltonian then takes the known central 
spin form H = J2i A l S e ■ h [9, [Io| . where S e is the sin- 
gle electron spin, and A4 = A\ip(fi)\ 2 the non-uniform 
coupling constant with ip(fi) the electronic ground state 
wave function at site fi . The reformulation of the central 
spin problem in terms of the KLH should be particularly 
useful for numerical evaluations. 

To continue with the general case, it is convenient 
to go to Fourier space and rewrite H n in Eq. ([2]) as 
H n = 4 J2qSq-Iq, where 1$ = J2j e~ % ^ ?i Ij is the Fourier 
transform of Ij , and similarly for S$. Since A is a small 
energy scale in our case, we may perform a Schrieffer- 
Wolff (SW) transformation in order to eliminate terms 
linear in A, and thereby integrate out the electronic de- 
grees of freedom. Keeping the lowest order terms in A 2 
of the SW transformation, we are left with an effective 
Hamiltonian H e ff = Hq — ^[S, [S, Ho]]. S is defined by 
H n + [S, Ho] = 0, which is solved as S = L ~ 1 H n where 
Lq is the Liouvillian. Let us define U — ^[S, [S, Hq]] 
which can be rewritten as U — \[Lq H n ,H n ], Us- 
ing an integral representation for L , one obtains U — 
— I J dte~ nt [H n (t), H n ], where 77 — > + ensures conver- 
gence. We next take the equilibrium expectation value 
over electronic degrees of freedom, denoted by (...). The 
only assumptions we make are (Sf) = (Sf) = 0, and 
translational invariance in the 2DEG. We then get 
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where Xap(q) = —i fo°dt e~ nt ({S~, S^_~]), and where 
summation over the spin components a, (3 — x,y,z is 
implied. If we also assume (Sf) = 0, then Xap(q) = 
$af3Xs(q), where Xs(q) is the electronic spin suscepti- 
bility in the static limit. We stress that Eq. ([3]) is 
rather general and requires only weak assumptions on 
Hq. In real space we have (U) = -jEf/^-^?^i 
where J"^ = — (A 2 /4n s )xa/3('r) is the effective exchange 
coupling. The nuclear spins Ip are therefore interact- 
ing with each other, this interaction being mediated by 
the conduction electrons. This is nothing but the stan- 
dard Ruderman-Kittel-Kasuya-Yosida (RKKY) interac- 
tion, which, as we shall see, can be substantially modified 
by e-e interactions compared to the free electron case. 

Let us first analyze the case of non-interacting elec- 
trons. In this case, Xs coincides with the usual density- 
density response (Lindhard) function xo . — — | - We first 
perform a mean field analysis. The Weiss mean field the- 
ory predicts a Curie temperature 
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where I is the nuclear spin value. In 2D, xo(q = 
0) = — N e = —m*/Tr, where N e — n e /Ep is the elec- 
tronic density of states, and m* is the effective elec- 
tron mass in a 2DEG (we set h = 1). For a 3D bulk 
metal with one conduction electron per nucleus, the ra- 
tio n e /n s r~j 1 and we recover the result in Eq. (jTJ) derived 
more than sixty years ago by Frohlich and Nabarro [13 ] . 
For a 2D metal, the Weiss mean field theory predicts 
k B T c = 1(1+1) A 2 IVlEp. For a 2D semiconductor, how- 
ever, the ratio n e /n s is much smaller than 1. With typical 
values for GaAs heterostructures, 7 = 3/2, A ~ 90 [ieV 
and a ~ 2A !)|. we estimate T c ~ 1 uK, which is very 
low. (For such low T c 's, ignoring nuclear dipole-dipole 
interactions from the start would not be valid.) How- 
ever, this estimate is just based on the simplest mean 
field theory and, moreover, does not include the effect of 
e-e interactions. 

We shall now go beyond above mean field approxima- 
tion. For this we assume that the ordering (if it takes 
place) leads to a ferromagnetic phase where the collective 
low-energy excitations are given by spin waves. Then, 
we define the Curie temperature T c as the temperature 
at which the magnetic order is destroyed by those spin 
waves. This procedure is equivalent to the Tyablikov de- 
coupling scheme (24|. The dispersion relation of the spin 
wave (or magnon) reads 
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I(J -J q )=I^a 2 (xs(q)-Xs(0)), (5) 
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where J q is the Fourier transform of J? . The magnetiza- 
tion m per site at finite T is m(T) = 7— j- ^-n g , where 
n q = [e^il k BTc _ 2)- 1 i s the magnon occupation number. 
The Curie temperature T c follows then from the vanish- 
ing of the magnetization, i.e. m(T c ) = 0, which, in the 
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continuum limit, becomes 
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For non-interacting electrons in 2D, Xs{q) — Xs(0) = 
for q < 2kp [22], where kp is the Fermi wave vec- 
tor. The spin wave analysis therefore predicts T c = 0, 
in agreement with a recent conjecture extending the 
Mermin-Wagner theorem for RKKY interactions in a 
non-interacting 2D system [251 ]. 

The study of thermodynamic quantities in interacting 
electron liquids especially in 2D has attracted quite some 
interest recently with the goal to find deviations from 
the standard Landau-Fermi liquid behavior, such as non- 
analytic dependences on the wave vector [TH, EE, G3, Ell- 
in particular, it was found that the static non- 
uniform spin susceptibility \s (<?) depends linearly on the 
wave vector q — \q\ for q — > in 2D (while it is q 2 in 3D). 
This non-analyticity arises from the long-range correla- 
tion between quasiparticles mediated by virtual particle- 
hole pairs. Since the integral in Eq. J6]) is dominated by 
the low q-behavior, one may replace u q by its low-g limit 
which turns out to be linear in q (see below) [23]. The 
integral in Eq. ((6]) can then be performed easily, allow- 
ing us to express T c in terms of the derivative of the spin 
sucseptibility, 
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For non- interacting electrons, Sxs(q) — at low q and 
we recover T c = 0, in accordance with the MWT. 

Let us include now e-e interactions. To calculate Xs(q), 
we start from the Bethe-Salpeter (BS) equation for the 
two-body scattering amplitude [22J. Solving the BS equa- 
tion formally, we can derive an exact and closed expres- 
sion for the spin susceptibility given by 
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where L = V Na 2 is the system length, (T ir )pp> (q) the ex- 
act irreducible electron-hole scattering amplitude in the 
spin channel (see [H), R p (q) = -2iG{p+ q/2)G(p- q/2) 
is the electron-hole bubble where G(p) is the exact prop- 
agator and p = (j>o,P) is the (D+l)-momentum with po 
the frequency. We have used a matrix notation in Eq. 
jH| where the indices run over p (R is a diagonal ma- 
trix). Unfortunately, T~ r cannot be calculated exactly 
and some drastic approximations are required. The ap- 
proximation we use consists in replacing the exact ir- 
reducible electron-hole scattering amplitude (r~ r )p,p> by 
an averaged value calculated with respect to all possible 
values of p and p' near the Fermi surface, therefore we 
assume {V~ r ) p ^ = T^ r (q) Vj),p' 

Let us now put qo = (and suppress the ^-argument 
from here on) and consider a ^-independent short-ranged 



(screened) interaction potential, yielding T~ r (q) = —U. 
This allows us to derive from Eq. (J5|) a simple formula 
for dxs/dq given by 



dxs_ , = dn{q) 
8q [q> dq 
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where H(q) = ~Ylp Rp{q) / L D . In the q — > limit, one can 
approximate the term II(<7) in the denominator of Eq. J9]) 
by Xo(0) = -N e . The resulting factor 1/(1 - UN e ) 2 in 
Eq. ([9]) can be interpreted as a type of random phase ap- 
proximation (RPA) for the electron-hole scattering am- 
plitude [l?]] ■ The corrections to the polarization bub- 
ble 11(g) (dominated by the first bubble correction to 
the self-energy) have been calculated in second order in 
perturbation theory (in U) at small q by Chubukov and 
Maslov The result of this perturbative approach 



is 6Tl(q) = U{q) - n(o) 
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r s ~ ~Um*/4n denotes the backscattering amplitude. 
When UN e <^ 1, we recover from Eq. ^ the known 
result 6 X s(q) = SU(q) 

Now we are ready to obtain an estimate for the Curie 
temperature T c . Replacing Xs(0) in 6xs(q) by its non- 
interacting limit xo(0), and assuming T s = 0(1) (this 
is an upper bound because T s is a small parameter con- 
trolling the perturbation theory), we obtain then from 
Eq. I© T c - 25 fiK for typical 2DEG parameters. This 
value of T c becomes further enhanced by a numerical fac- 
tor (e.g. of order 5 for r s ~ 8 [22]) if one uses an effective 
renormalized value for the spin susceptibility xs = Xs(0) 
instead of Xo(0). Though T c is still rather small, it is now 
finite, confirming our arguments related to the Mermin- 
Wagner theorem that e-e interactions increase the Curie 
temperature. When UN e is no longer negligible com- 
pared to 1, T c is even further enhanced by an additional 
numerical factor 1/(1 - UN e ) 2 (see Eq. ©). Close to 
the ferromagnetic Stoner instability of the electron sys- 
tem, reached when UN e ~ 1, the Curie temperature T c 
for the nuclear system is dramatically enhanced as could 
have been anticipated. 

In the preceeding paragraphs, we replaced T~ r (q) by a 
q— independent constant operator. One can use instead 
another approximation called the local field factor ap- 
proximation (LFFA) . The idea of the LFFA is to replace 
the average electrostatic potential by a local field poten- 
tial seen by an electron with spin a (see [12] for a re- 
view). In this scheme (r~(<?)) P j/ « — V(q)G- (q), where 
G_ (q) is a local field factor and V(q) = 2ire 2 / nq the bare 
unscreened Coulomb interaction (k is the dielectric con- 
stant). Within this approximation scheme the static spin 
susceptibility Xs becomes 
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Determining precisely G-(q) for all q is still an open 
issue. However, the asymptotic regimes are quite well 
established nowadays [23]. A semi-phenomenological in- 
terpolation formula based on the original Hubbard local 
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field factor [28J] and modified in such a way that the com- 
pressibility sum rule is exactly satisfied reads [H, 



G- (q) ~ 5o 
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where go is related to the probability of finding two elec- 
trons (of opposite spins) at the same position in the elec- 
tron liquid, (gfiB)~ 2 XP ls the Pauli susceptibility and 
fiB the Bohr magneton. For non- interacting electrons 
Xp/xs — 1- An approximate form for go giving good 
agreement with quantum Monte Carlo (QMC) calcula- 
tions has been proposed recently by Gori-Giorgi et al. 
[Ij: g {r s )~(l + Ar s +Br 2 s +Cr 3 s )e- Dr °/2. In a 2DEG, 
r s = 1/ \jim e a* B where a* B — n/m*e 2 is the effective Bohr 
radius. The parameters A = 0.088, B = 0.258, C = 
0.00037, D — 1.46 are fit ting parameters reproducing 
QMC results for the 2DEG (if. From Eqs. ([7]) and (fi"Q|) . 
one can easily determine T c within the LLFA scheme to 
be given by 
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where a = (1 — xp/xs)^ 1 and V(a) is the Coulomb po- 
tential evaluated at the interatomic distance a. The en- 
ergy scale (a — l) 2 goV(a) can be interpreted as a renor- 
malized screened potential due to collective interaction 
effects that are incorporated in the LFFA. The ratio 
A/(a — l) 2 goV(a) can be regarded as the small param- 
eter of our theory. Quite remarkably, the LFFA pre- 
dicts an exponential enhancement of T c with increasing 
interaction parameter r s . For a value of r s ~ 5, this 



theory already predicts a large T c ~ 25 mK, a tem- 
perature which is routinely achieved nowadays. Obvi- 
ously, for some value of r s , the dimensionless parameter 
A/(a — l) 2 goV(a) exceeds unity. The truncation of the 
Schrieffer- Wolff transformation at lowest order becomes 
unjustified and feedback effects between the electron gas 
and the nuclear spins, not incorporated in our theory, be- 
come important. Nevertheless for relatively large values 
of r s < 6, the condition A <§; (a — l) 2 goV(a) is satisfied. 

Although the spin wave analysis may overestimate T c , 
the trend in all the approximation schemes we used is 
that e-e interactions increase dramatically the Curie tem- 
perature, possibly into the mK range for large r s (there- 
fore three orders of magnitude larger than Edi P which 
justifies our starting Hamiltonian) . We note that the 
non-perturbative LFFA theory predicts higher T c 's than 
the perturbative calculation in the short-ranged inter- 
action. Finally, below T c , the nuclear spins within the 
2DEG polarize and generate an effective magnetic field 
of order of a few Tesla. This will create a small Zeeman 
splitting [1H in the 2DEG which should be detectable 
with e.g. optical or transport methods. 

In summary, we have analyzed the Curie tempera- 
ture T c of nuclear spins in an interacting 2DEG using 
a mean field and a spin wave analysis. We have shown 
that electron-electron interactions considerably enhance 
the temperature for a ferromagnetic phase transition in 
the nuclear system, with T c in the milli-Kelvin range for 
2DEGs with r s ~ 5 - 10. We thank B. Coish, L. Glaz- 
man, L. Kouwenhoven, and A. Yacoby for useful discus- 
sions. This work is supported by the Swiss NSF, NCCR 
Nanoscience, ONR, and JST ICORP. 
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